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size,  and  mass  of  the  particles.   The  work  treats  the  mass  as  a  first  order 
perturbation,  in  the  particle  momentum  equation.   The  unsteady  flow  studied  is 
a  region  in  the  wake  of  a  cylinder  where  charged  water  droplets  are  injected 
into  an  air  flow.   From  this  work  the  size  of  the  charged  particles  is  calcu- 
lated and  a  correspondence  is  found  between  the  mean  particle  size  injected 
and  the  observed  frequency  response  of  the  anemometer  tested.   The  EGD  spectral 
anemometer  data  are  compared  to  hot-wire  data  taken  at  the  same  location;  some 
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frequency  response  as  a  function  of  charged  particle  size  are  given. 
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I.   INTRODUCTION 

Electrof luiddynami  c  devices  operate  on  the  principle  that  a 
nonconducting  fluid  can  interact  with  an  electric  field  when  unipolar 
charges  are  introduced  into  the  fluid.   If  this  interaction  does  not 
disturb  the  fluid  appreciably  and  if  the  charged  particles  are  linked 
intimately  with  the  flow,  then  the  arrival  of  these  particles  at  a 
suitable  collector  brings  information  about  characteristics  of  the  flow. 
This  work  deals  exclusively  with  gases  and  the  term  electrogasdynamic 
(EGD)  anemometer  will  be  used  to  denote  the  device  under  consideration. 

The  feasibility  of  the  anemometer  hinges  on  the  proper  choice  of 
charged  particles.  These  must  couple  the  gas  to  the  electric  field 
through  individual  collisions  if  the  particles  are  small  enough  or 
through  fluid  dynamic  drag  and  pressure  forces,  if  they  are  large  enough. 
Molecular  ions  are  too  mobile  to  be  useful  ;  therefore,  charged  particles 
of  size  much  larger  than  the  gas  molecules  must  be  introduced.  The 
ultimate  size  of  these  particles  is  restricted  by  their  ability  to 
follow  the  turbulent  fluctuations  and  by  the  fact  that,  depending  on  the 
method  of  charging,  a  large  particle  may  have  a  sizeable  mobility  (as 
will  be  seen  in  Section  III).   In  this  paper,  the  effect  that  particle 
size  and  charge  have  on  the  coupling  between  the  fluid  and  the  electric 
field  will  be  discussed. 

It  may  be  mentioned  in  passing  that  the  definition  of  particle  slip 
in  unsteady  flow  is  not  only  crucial  for  the  anemometer  but  also  for  the 

EGD-power  generator.   Since  turbulence  is  a  means  of  enhancing  the  break - 

2  3  1  h 

down  strength  in  a  gas  '   and  since  generators  '  rely  on  the  mixing  of 

fluids  both  to  increase  the  pressure  drop  in  the  channel  and  to  insulate, 

there  must  exist  regions  of  unsteadiness  in  the  flow.   Hence,  a  small 


slip  parameter  for  steady  flow  may  misrepresent  the  ability  of  a  charged 
particle  to  interact  with  the  flow. 
II.   EGD  -  Anemometer 

V7hen  a  charged  particle  is  tightly  coupled  to  the  flow,  it  is 
essentially  a  flow  "tracer."  To  pick  up  such  a  tracer  one  simply  places 
a  metallic  probe  which  provides  the  neutralizing  current.  The  injected 
charges  need  a  certain  time  to  equilibrate  with  the  flow  but  there- 
after they  drain  no  energy  from  the  flow.     The  particles 
are  neutralized  as  they  arrive  with  the  surrounding  gas  and  are  not 
forced  to  the  probe  by  an  external  field  hence,  in  principle,  the  probe  does  not 
distrub  the  nature  of  the  flow. 

Because  the  collection  area  is  not  necessarily  the  probe  area  and  be- 
cause of  an  uncertainty  in  the  size  and  charge  of  the  particles, 
it  is  expedient  to  measure  the  spectral  characteristics  of  the  flow 
independently  with  a  hot-wire  anemometer  and  then  to  compare  the  results 
with  those  of  the  EGD -anemometer .  This  permits  the  optimization  of  the 
operating  parameters  of  the  EGD -anemometer  and  allows  the  inference  of 
charge  and  size  distributions  for  the  charged  particles. 

The  hot-wire  anemometer,  although  a  heat  transfer  measurement,  has 
emerged  as  perhaps  the  most  useful  device  for  measurements  in  unsteady 
flows.   Needless  to  say,  advances  in  fluid  mechanics  hinge  on  reliable 
measurements,  and  hot-wire  data  are  an  accepted  standard  for  these 
measurements.  The  EGD  device  can  complement  the  hot-wire  anemometer. 

An  EGD -anemometer  consists  of  the  following:   an  injector  or  source 
for  the  charged  particles,  a  conversion  region  where  the  particles 
equilibrate  with  the  flow,  and  a  collector  for  the  charges.   The  collector 
unit  corresponds  to  the  anemometer  probe.   The  EGD-probe  is  very  rugged, 
allowing  the  probing  of  intense  turbulent  fields.   It  is  simple  and 


inexpensive,  and  therefore  can  be  easily  manufactured.   Since  very  small 
charged  particles  can  react  to  high  frequency  oscillations,  the  frequency 
response  of  the  EGD- probe  can  be  as  high  as  the  best  hot-wire  (at  no 
compromise  in  mechanical  strength) .   In  addition,  signals  from  the  EGD- 
probe  are  electrical  and  need  little  or  no  further  amplification.  Effects 
corresponding  to  the  thermal  inertia  of  hot-wires  are  not  present.   There 
is  no  need  for  a  linearizer  since  the  EGD  measurement  is  already  linear,  i.e., 
a  direct  measurement.  Absent  are  the  complications  that  free  convection 
currents  bring  to  the  hot-wire  zero  setting.   In  the  case  of  two- phase  flows, 
the  EGD-probe  needs  no  separate  calibration  and  will  withstand  the  impact 
of  charge  droplets  or  solid  particles  whereas  the  fragile  hot-wire  will 
readily  break  under  these  conditions.   In  the  case  of  two  mixing  jets, 
there  will  invariably  be  a  temperature  difference  and  this,  however  slight, 
will  affect  the  hot-wire  since  it  is  sensitive  to  temperature  differences. 
The  geometry  of  the  EGD  probe  is  more  adaptable  than  that  of  the  hot-wire 
to  unusual  probing  situations. 

The  main  disadvantage  of  the  EGD  scheme  for  anemometry  is  the  require- 
ment for  a  source  of  particles  of  optimum  size  and  charge.   If  the  particles 
are  created  through  a  corona  discharge  then  there  is  a  need  for  a  high 
voltage  supply  (with  the  associated  safety  precautions).   Also,  the 
injector/ collector  unit  must  have  a  separation  of  the  order  of  millimeters 
and  must  be  traversable  all  of  which  pose  several  design  problems.   However, 
all  these  problems  can  be  overcome  and  the  use  of  an  EGD-generator  scheme 
as  a  spectral  anemometer  looks  quite  feasible. 

It  should  be  pointed  out  that  the  use  of  charges  in  a  flow  for 
anemometry  is  not  entirely  a  new  idea.  Hinze  describes  electric  discharge 
anemometers  which  consist  typically  of  a  glow  or  corona  discharge  in  the 
flowing  medium.  More  recently,  Nygaard  has  worked  on  a  type  of  corona- 
discharge  anemometer.  A  mass  flow  meter  device  based  on  an  ion- displacement 

7 
concept  of  Professor  Durbin  is  being  marketed  by  Thermo- Systems,  Inc. 


None  of  the  reported  schemes,  however,  have  concentrated  on  the  spectral 
response  of  the  instrument  and,  in  particular,  in  matching  the  charged 
particles  to  the  desired  value  of  the  slip  at  a  given  frequency  range. 
III.   Concepts  from  Unsteady  Flow 

The  ability  of  a  charged  particle  to  couple  its  surrounding  gas  to 
the  electric  field  is  described  by  its  slip.  The  drift  velocity  with 
which  a  charged  particle  moves  in  an  electric  field  is  given  by  the 
product  of  its  mobility  (a.)  times  the  electric  field  (E) 

vD  =  u.E  (1) 

and  the  slip  (S)  for  steady  flow  is  defined  as 

S  "  T 

V  =  fluid  velocity 

Q 

Here  one  wishes  to  make  S  as  small  as  possible. 

So,  knowing  the  electric  field  and  the  fluid  velocity,  the  slip  can 
be  calculated  provided  that  the  mobility  is  also  know.   Unfortunately, 
deducing  the  mobility  is  not  straightforward  since  it  depends  on  the  size 
of  the  particle  and  on  how  it  was  charged  .   If  one  plots  the  drift 
velocity  defined  in  Equation  1  for  the  maximum  field  (breakdown)  as  a 
function  of  the  radius  for   spherical  water  droplets  in  air  at  standard 
conditions,  one  gets  the  curves  shown  in  Figure  1.  Singly  charged 
particles  have  a  decreasing  [iK   -velocity  with  radius  because  the  drag  is 
increasing  while  the  electrical  force  remains  constant.   The  case  for 
particles  charged  to  the  Rayleight  limit  represents  an  upper  bound  in  the 
charging  scheme  because  more  charges  exceed  the  stability  of  the  droplet 
at  any  given  radius.   Corona  field  charging  is  an  intermediate  case 
between  the  other  two  modes  of  charging.  ■  Since  the  speed  of  sound  is 
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FIGURE    I. 


MAXIMUM   DRIFT  DUE  TO  ELECTRIC  FIELD   FOR 
WATER   DROPLETS    IN  AIR  AT   I -ATM   AND   20°  C 


representative  of  flow  velocities  of  practical  interest,  Figure  1  gives 
information  about  the  range  of  radii  and  the  charging  scheme  that  would 
be  required  for  a  given  value  of  the  slip. 

When  the  flow  is  unsteady,  the  question  arises  as  to  what  happens 
to  the  conventional  slip  parameter.   Clearly,  inertial  effects  of  the 
particles  will  come  into  play.   In  Appendix  A,  the  momentum  equation  for 
a  particle  is  treated.  There  a  one -dimensional,  spectral  analysis  is 
given  which  shows  that  the  concept  of  the  slip -parameter  may  be  extended 
to  unsteady  flow.  Written  in  a  spectral  form,  the  maximum  value  of  the 
slip  becomes 

max   V(o)')     v 


B(»)|„  =  yfa)  +UJT"  (2) 


where 


uj  =  frequency 

V(ou)  =  RMS  spectral  intensity 

t  =  e -folding  time  for  the  particle  velocity  to  approach 
steady  state  value 

t   (the  equilibration  time)  is  a  function  of  the  particle  radius  and 

of  the  viscosity  of  the  fluid.   If  turbulence  or  any  unsteady  flow  is  to 

have  an  effect  on  the  charges,  then  the  slip  parameter  must  be  small  for 

all  the  dominant  frequencies  (tu  )  of  the  flow,  i.e., 

HE 

+  UJ,   T   «  1.0 


vkT  wk  v 

k 
Actually,  it  is  only  necessary  to  look  at  the  highest  frequency  of 

interest  since  for  all  smaller  values  of  uo,  the  inequality  above  will  be 

satisfied  provided  V(cu  )  retains  nearly  the  same  magnitude. 

Figure  2  shows  two  curves  bracketing  a  region  in  which  the 

-3       -2 

contribution  of  urr  to  Equation  2  is  between  10   and  10   .   This  range 

represents  a  criterion  of  small  slip.   The  contribution  of  the  first 
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FIGURE     2.        OPTIMUM       FREQUENCY    RESPONSE 

AS     A     FUNCTION     OF    THE     RADIUS 


term  in  the  right-hand  side  of  Equation  2  may  or  may  not  be  as  low  as 

-2 

10   hut  one  would  endeavor  to  pick  a  small  enough  value  by  manipulating 

with  the  charge  of  the  particles  (see  Figure  l) . 

In  Figures  1  and  2,  water  droplets  have  been  described. 
For  less  dense  particles  or  for  hollow,  spherical  solids  the  curves  move 
to  the  right  (i.e.,  towards  larger  radii)  which  is  a  desirable  effect 
from  practical  considerations.   It  is  conceivable  that  one  size  of 
particle  will  be  optimum  for  many  applications. 

The  measurement  of  particle  size,  or  more  accurately  stated  of 

particle  size  distribution,  is  very  difficult  in  the  range  of  interest 

9  10 

of  EGD-anemometry  '   .  The  water  droplets  can  only  be  measured  in  situ 

and  direct  techniques  such  as  laser  holography   can  presently  not 
measure  below  50  x  10  m*.   In  the  course  of  this  work  an  indirect 
method  of  measurement  was  developed  which  consists  of  comparing 
the  spectrum  of  the  current  collected  to  that  of  the  hot-wire 
spectrum.  For  the  proper  particle  size  and  charge  distribution,  the 
agreement  between  the  two  measurements  is  high  when  the  flow  has 
not  been  significantly  disturbed  by  the  particles. 
The  current  collected  may  be  written  as 

00 

Ig(m)  =  So   vpKR)q(R)n(R)dR  (3) 

where  R  =  radius  of  the  particle. 

The  total  number  of  particles  per  unit  volume,  n  ,  is 


00 


n(R)dR  =  n 
Jo   ^  '      1 


o         p 

and  q(R)  =  charge  of  the  particle 

v  (uu,R)  =  velocity  of  the  particle 


Dividing  Equation  3  by  v(uu),  the  flow  velocity,  which  is  independent 

of  the  particle  radius,  one  gets 

I  (cu)    «>  v  (uu,R) 

-4-t-  =  f  Vl    \        q(R)n(R)dR  (k) 

It  is  shown  in  Appendix  B  that  the  ratio  v  (uj,R)/v(oj)  may  be  approximated 
by  a  Heaviside-type  of  function, 

v  (w,R) 

P/  .   =1  t  <  1/uu 

v(cju)  V    ' 


=  0  T   >  1/ 


CO 


V 


This  defines  R'(cd),  the  radius  at  which  t  equals  cju   and  Equation  k 


becomes 


I  (cu)    R'(co) 
-vW^Jo     *<*)n(R)dR 

where  for  spherical  particles 

(j,  =  viscosity  of  fluid 
p  =  particle  density 
Since  one  can  obtain  I  (tu)  and  v(uu)  experimentally,  and  since  one 

o 

can  surmise  what  q(R)  must  be,  one  can  deduce  the  distribtuion  n(R)  by 

solving  the  implicit  equation  above  by  successive  approximations.   q(R) 

is  given  by  the  radius  of  the  particle  and  the  method  of  charging,  as  can 

be  seen  from  Figure  1.   If,  however,  the  spread  of  particle  size  is 

narrow  enough,  then  an  average  q  per  particle  may  be  used  as  follows: 

I  (cju)      R'  (cju) 
W  =  *  Jo    "(R)dR 

or  (5) 

g  ,  '      P  =  ±-   J     n(R)dR 
v(oj)      n  Jo     v  ' 

P 
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Therefore,  the  ratio  of  the  "adjusted"  value  of  the  current  to  the 
fluid  velocity  represents  the  integral  of  the  particle  size  distribution. 
When  the  radius  R'(ou),  equivalent  to  a  particular  frequency  in  question, 
is  larger  than  that  of  the  particles,  the  ratio  of  Equation  5  will  be  a 
constant.  When  it  lies  within  the  range  of  the  particle  size  distribt- 
tion,  however,  a  decreasing  value  for  this  ratio  will  be  seen  as  the 
frequency  increases.  At  a  high  enough  frequency,  corresponding  to  radii 
smaller  than  the  particles,  the  ratio  will  go  to  zero. 

The  hardware  and  the  measurements  are  treated  next  before  Equation  5 
is  further  discussed. 
IV.   Experimental  Set -Up 

In  the  EGD  apparatus,  water  droplets  are  introduced  into  an  air- 
streams  as  depicted  in  Figure  3-  The  injector  consists  of  an  aerosol 

flow  in  conjunction  with  an  ionizer.  Saturated  steam  is  injected  through 

12 
a  corona  discharge   into  the  turbulent  wake  of  a  cylinder.  Electrons 

from  the  needle  drift  toward  the  attractor  and  attach  themselves 

to  the  then-growing  water  droplets  forming  negative  charged  particles . 

This  arrangement  with  the  steam  as  the  source  of  particles  permits  some 

variation  of  particle  size  distribution  by  mere  manipulation  of  the  steam 

pressure  and  temperature. 

The  anemometer  is  tested  in  a  plexiglass  wind-tunnel.  The  inside 
dimensions  are  five  cm  by  ten  cm  in  cross  section  and  25  cm  in  length. 
The  test  section  contains  the  steam  injector  unit.  The  aerosol  injector 
nozzle  consists  of  a  stainless  steel  ring  2.7  mm  in  diameter.  This 
arrangement  is  shown  in  Figure  K. 

The  steel  ring  also  serves  as  the  corona  ring  positive  electrode 
(attractor)  of  the  aerosol  ionizer.   The  negative  electrode  consists  of 
a  stainless  steel  needle.  Typically  at  five  millimeters  downstream  from 
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FIGURE     4.         AEROSOL      INJECTOR 
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the  cylinder  nozzle,  a  third  electrode,  in  the  form  of  a  slender  rod 
1.5  mm  diameter,  collects  the  current.  This  separation  distance  is  far 
enough  downstream  for  equilibration  between  the  particles  and  the  flow  to 
take  place  and  close  enough  for  substantial  currents  to  be  collected. 

The  air  flow  is  fed  through  a  temperature  regulation  stage  for  the 
control  of  the  relative  humidity  since  it  was  found  that  water  droplets 
were  somewhat  affected  by  the  humidity  of  the  air  during  mixing  outside 
of  the  nozzle.  The  humidity  of  the  air  was  held  at  32$  because  results 
were  most  reproducible  at  this  value.  (Humidity  control  would  not  be 
necessary  with  solid  particles.)  The  water-air  mixture  flows  out  of  the 
test  section  into  the  atmosphere. 

The  output  from  either  the  hot-wire  or  the  EGD -collector  probe  is 
fed  into  a  black  box  which  houses  the  circuitry  shown  in  Figure  5«  The 
purpose  of  this  black  box  is  primarily  for  the  protection  of  instrumenta- 
tion downstream  of  the  probe.  The  spectral  data  are  obtained  from  a 
frequency  analyzer.   It  is  a  General  Radio  Company  Type  1920  Real  Time 
Analyzer.   It  consists  of  a  multifilter  band-pass  section  and  a  multi- 
channel RMS  Detector.  The  signal  from  the  preamplifier  is  fed  into  the 
multifilter,  separated  according  to  frequency,  then  fed  into  the  RMS 
detector  where  the  data  are  converted  to  binary  digital  form  and  the  RMS 
level  computed  by  true  linear  integration. 

Hot-wire  units  of  the  constant -temperature  type  with  frequency 
responses  over  50  KHz  were  utilized.   The  hot-wire  readings  were  used  as 
a  standard  of  comparison  for  the  EGD  measurements  but,  because  the  hot- 
wire calibration  is  not  appropriate  for  two-phase  flows,  air  was 
substituted  for  the  steam  in  the  cylinder.  The  air  flow  was  matched  so 
as  to  give  the  same  momentum  as  the  steam.   The  absence  of  droplets  also 
helped  prolong  the  life  of  the  hot  wires  which  was  very  marginal  to  begin 
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with  in  the  intense  vorticity  of  the  wake  of  the  cylinder.   The 
resulting  hot-wire  profile,  after  the  temperatures  of  the  two 
jets  were  equalized,  is  shown  in  Figure  6  for  both  the  main  air  and 
nozzle  air.   There  are  two  dominant  frequencies  apparent  on  this 
figure,  namely,  1.2  KHz  and  lOKHz.   The  first  one  has  the  value  of 
twice  the  Strouhal  frequency  for  the  cylinder  and  is  appropriately 
so  because  measurements  were  taken  along  the  centerline  of  the 
wake  where  the  vortices  shed-off  by  the  cylinder  on  alternating  sides 
add  up.   The  Reynolds  number  based  on  the  cylinder  diameter  for  this 
flow  (56  m/sec)  is  8  xllO  so  that  the  boundary  layer  on  the  cylinder 
is  laminar  before  separation.   The  second  important  frequency  in  Figure  6 
is  associated  with  the  nozzle  for  the  steam-- the  inner  jet.   The  10  KHz 
peak  appeared  whether  or  not  there  was  air  issuing  from  this  nozzle.   The 
exit  diameter  for  this  nozzle  is  2  mm  which  is  some  ten  times  smaller  than 
the  cylinder  diameter  and  appropriately  the  frequency  is  nearly  ten  times 
higher.   It  was  observed  that  10  KHz  is  a  high  enough  frequency  so  that 
the  dissipation  of  the  vortices  is  high;  at  a  station  10  mm  away  from 
the  nozzle  exit  the  10  KHz  peak  disappears  whereas  the  1.2  KHz  peak 
remains  nearly  as  intense. 
V.   Results  and  Discussion 

The  operation  of  the  EGD- anemometer  was  investigated  for  frequency 
response  and  for  current  convection  efficiency.   The  high  frequency  response 
was  most  noticeable  when  the  steam  was  very  close  to  the  saturated 
condition  before  expansion  (charges  of  1  C  noticeably  affected  the  results). 
Because  of  this,  it  was  decided  to  run  the  steam  at  two  stagnation 

pressures  to  achieve  different  size  of  particles.   Previous  work  in  our 

13 

laboratory   indicated  that  for  the  aerosol  injector  used  and  for  an  air 

flow  of  56  m/sec,  smaller  droplets  appeared  at  steam  reservoir  pressures 

below  1.6  atm. 
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The  method  of  comparison  between  the  EGD  and  hot-wire  data  is  based 

on  Equation  5.      T   (uu)   is  defined  as  an  RMS  parameter,   namely, 
6 

I     /An  q 

where 

I   =  spectral  probe  current  (RMS) 

A  =  collector  area 
n  =  total  number  of  charge  particles  per  unit  volume 

q  =  average  charge  per  particle 
v  =  flow  velocity  (RMS) 

Letting  u  be  a  reference  voltage  (one  millivolt  in  this  case),  one 

obtains  for  the  value  of  T  in  dbs , 

g 

I  /An  q  v 

20  log..  _[T  (ou  )]  =  20  logn^  -ffl 2-  -  20  log,n  — 

10  g  n         °10    u  10  u 


(6) 


T  (ou  )  T^u)  )  T,  (u)  ) 

gv  n7  Ev  n'  hv  n' 


where 


T  (ou  )  =  measurement  of  the  EDG-probe  in  dbs 

T,  (u>  )  =  measurement  of  the  hot-wire  in  dbs 
h  n/ 

The  points  shown  in  Figure  7  through  10,  have  the  original  EGD -probe 

data  and  a  curve  corresponding  to  the  difference  between  those  data  and 

the  hot-wire  values  of  Figure  6.   Twice  the  Strouhal  frequency  was  chosen 

as  a  reference  point  in  this  comparison  because  it  is  the  most  intense 

point  and  because  frequencies  higher  than  1.2  KHz  are  most  sensitive 

to  changes  in  the  steam  parameters. 
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Figures  7  and  8  show  typical  results  for  a  steam  reservoir  pressure 

of  1.85  atm  but  with  differing  temperatures  and  corona  voltages.  At 

I.85  atm,  the  corresponding  saturation  temperature  is  118  C  whereas  the 

measured  values  were  116  C  and  117  C  for  the  data  on  Figures  7  and  8 

respectively.  Figures  9  and  10  show  typical  results  at  1.55  atm.  Here 

the  saturation  temperature  is  113  C  so  it  is  clear  that  some  super - 

saturation  was  achieved  at  both  operating  pressures  before  the  expansion. 

If  the  curve  for  T  were  a  horizontal  line,  then  the  two  sets  of  measure- 

g 

ments  would  only  differ  by  a  constant  and  the  EGD  anemometer  would 
essentially  reproduce  the  hot-wire  anemometer  measurements.  An  examina- 
tion of  Figures  7  through  10  shows  that  the  present  EGD-probe  can  measure 
intense  turbulent  fields  with  some  degree  of  success.  Figure  8  shows  a 
correspondence  of  within  +  2dbs  up  to  630  Hz  and  Figure  9  shows  a 
correspondence  within  +  6dbs  up  to  20  KHz.  The  EGD  -probe's  response  can 
be  improved  by  operating  with  droplets  of  considerably  smaller  size 
distribution  as  will  be  discussed  below. 

The  technique  outlined  earlier  in  this  paper  is  entirely  feasible 
for  the  measurement  of  particle  size  distribution.   It  is  evident  from 
Figures  7  and  8  that  the  data  at  I.85  atm  are  well  represented  by  the 
assumptions  stated  in  Section  III.  There  it  was  said  that  when  the 
particle  size  distribution  is  such  that  it  overlaps  into  the  region  of 
interest  in  the  frequency  scale  (from  Equation  5)  then  the  amplitude  of 
the  signal  will  diminish  in  a  monotonic  fashion.  The  data  in  Figure  7 

may  be  interpreted  as  indicating  a  distribution  with  center  radius  of 

-5  -5 

2  x  10   m  and  with  a  half -width  of  about  1.5  x  10   m  (assuming  a 

Gaussian  distribution  for  the  charged  droplets  and  a  corresponding 

average  charge  per  droplet).   The  data  in  this  figure  have  a  correspondence 

with  the  hot-wire  only  up  to  about  80  Hz,  which  means  that  R'(uj)  is 
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2.9  x  10   m.  R'(uu)  is  essentially  the  higher  value  of  the  radius  at 

the  half -width.  The  data  from  Figure  8  has  a  correspondence  with  the 
hot-wire  up  to  630  Hz  and  this  yields  and  R'(uo)  of  1.2  x  10"^  m.  The 
data  at  1.55  atm,  while  having  a  better  overall  frequency  response,  do 
not  show  a  monotonic  decrease  with  increasing  frequency.  This  implies 
the  presence  of  more  than  one  Gaussian  distribution  for  the  particles 
(the  possibility  of  droplet  fragmentation  in  the  wake  of  the  cylinder  has 

TO 

been  proposed  ).   It  is  evident  that  conditions  outside  of  the  nozzle 
strongly  affect  the  droplet  size  and  the  statistical  sample  of  data 
available  at  the  present  is  too  sparse  to  allow  further  conclusions.   It 
can  be  seen  from  Figure  2,  however,  that  a  more  favorable  particle  size 
distribution  (i.e.,  centered  about  10   m)  would  enhance  the  EGD-probe's 
response  to  match  that  of  the  hot-wire  up  to  20  KHz. 

Turning  now  to  other  operating  characteristics  of  the  anemometer, 
Figure  11  shows  the  effect  of  the  density  of  the  steam  on  the  current 
convection  efficiency  (I  /i™)  for  the  two-phase  fluid.   In  addition  to 
the  sensitivity  of  this  efficiency  to  small  changes  in  voltage,  the  effect 
of  an  increase  in  density  is  to  require  a  higher  voltage  for  the  same 

efficiency.   This  shifting  of  the  curves  to  the  right  is  similar  to  that 

Ik 

reported  by  Whitby   for  the  effect  of  pressure  (temperature  was  held 

constant)  in  a  single-phase  fluid  (air).  Steam  saturation  effects  are 
not  apparent  when  pressure  and  temperature  data  are  converted  to  density 
but  they  play  an  important  role . 

Figure  12  shows  that  a  definite  optimum  voltage  exists  for  operation 
with  a  good  high  frequency  response.   After  a  "triggering"  voltage  is 
attained  the  collector  current  increases  linearly  until  finally  a  plateau 
is  reached.   The  graph  shows  that  the  optimum  current  efficiency  does  not 
occur  at  optimum  voltage  although  some  preliminary  observations  indicated 
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that  acceptable  results  occurred  at  a  current  convection  efficiency  of 
0.2  to  0.5  with  higher  collector  currents.   Figure  13  show:;  the  effect 
of  efficiency  on  the  high  frequency  response  for  three  different 
temperatures  with  other  parameters  held  constant. 
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Appendix  A 
One -Dimensional,  Spectral  Analysis 
For  the  case  where  the  energy  exchange  between  the  flow  and  the 
electric  field  is  slight,  only  the  particle  momentum  equation  is 
considered.  The  one-dimensional  case  depicted  in  Figure  lU  is  a  useful 
simplification --it  represents  an  extreme  case  where  the  entire  fluid 
velocity  is  opposite  to  the  electric  field.   In  turbulent  flows,  this 
would  be  true  only  for  part  of  the  time.  But  one  is  interested  here  on 
an  estimate  of  the  maximum  slip  and  hence  on  a  description  of  the 
maximum  interaction.   The  particle  analyzed  is  assumed  to  be 
large  enough  to  experience  fluid  dynamic  drag  forces  (i.e.,  its  diameter 
is  greater  than  the  mean  free  path  of  the  gas  molecules)  but  small 
enough  so  that  pressure  forces  may  be  neglected.  Stoke' s  drag  is  assumed 
applicable  even  though  this  overestimates  the  actual  drag  for  Reynolds 
numbers  (based  on  the  particle  diameter)  which  exceed  one.  In  a  restricted 

sense  the  particles  size  will  be  allowed  to  range  to  below  the  mean 

15 
free  path  of  the  gas 

The  free  stream  velocity  "v"  is  divided  in  to  spectral  components  at 

preselected  frequencies  co  ,  cu  ,  (i)  ,  . . .  a>  as  follows:   Let  the  velocity 

at  a  given  location  be 


v(t)  =  £V(w  )f(uj  ,t) 


where 

V(ou  )  denotes  the  amplitude 

and  f(w  ,t)  is  a  periodic  function  with 


n- 


|f(u>  ,t)|  =  1.0 
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Now,  at  a  given  frequency  uu  ,  the  particle  momentum  equation  can  be 

written  as 

dv  (uu  ,t) 
mp  ~^t1 =  6TTB^j:v(ci)n,t)  -  vp((i)n,t)]  -  qE  (Al) 

Stoke  's  Drag        Electric  Force 
m  =  mass  of  the  particle 

R  =  radius  of  the  particle  (spherical  shape) 
H  =  viscosity  of  the  fluid 

The  electrical  force  is  assumed  steady  here.  There  may  be  space  charge 
fluctuations  which  influence  "E"  locally  but  these  are  assumed  to  be 
outside  of  the  frequency  range  of  interest. 

The  particle  velocity  can  be  similarly  written  as 


v  =  E  V  (u)  )g(cu  ,t) 
p    n  p  n    n' 

where  V  is  the  amplitude  and  again  g(uu  ,t)  is  periodic  with 


|g(o)  ,t)|  =  1.0  so  that  Equation  Al  becomes 


dg(u)  ,t)   m 
Vp("n)  -dt =  ^YK)fK^    ~   V»n)«(»n.*)]  -  <& 

(A2) 

where 

m 

\ s  zira  (A3) 

which  is  the  equilibration  time  for  a  constant  velocity  and  no  electric 
fields. 

Now,  if  the  particle  mass  (m  )  is  vanishingly  small,  one  gets  the 
usual  definition  for  the  mobility  from  Equation  Al,  namely 


vD  =  v  -  vp  =  p-E  (Ak) 
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where 


v_  =  drift  velocity 
[i   =  particle  mobility  =  £-§ — 

Before  considering  the  inertia  term  as  a  perturbation,  Equation  A2 
will  be  nondimensionalized.   Let 

t  =  tco 
n     n 

so  that  for  the  uu     particle,   in  simplified  notation,   one  gets 


dg         m 

mVo)     -rt£-  =  r^Cv     -  v     ]    -  qE 
p  p  n  dt         t         n         pn 


n         v 


where 


or 


v     =  v(u)    ,t)   and  v       =  v  (u>    ,t) 
n  v   n'  pn         p     n'    ' 


.m  V  U)  .  dg  m 

(JLJLR) -*°  =  -^[v     -  v    ]    -  1  (A6) 

\      qE     /dt  t   qE       n         pn  v      ' 


_> 


0(O      0(1)  0(1) 

From  the  above  equation,   one  can  see  the  perturbation  parameter 
explicitly,   namely, 


m  V  en 

qE 


e    =     PP" 


Now  let 


Co)  +ea  CD  +0(c2' 


g     =  g   v    ;    +  eg    v    y   +  0(e    ) 

Substituting  into  Equation  A6,    collecting  terms,   and  neglecting  those  of 
0(e    ) ,  one  obtains 


r-^r  [v     -  V     g   (o)]    -1  =  0  (A7) 

t   qE        n         pnton  v      ' 
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and 


m 


.  (°)]  =  .  -je_rv  k  (1)i 

n  v    * 

From  the  second  of  the  two  equations  (Equation  A8) 


(A8) 


g 


(1)     _ 


t   qE     d| 


(o) 


n 


m  V 
P  Pn 


n 


Because  E  is   assumed  steady,   Equation  A7  replaces  Equation  A5, 
which  defined  the  mobility.      Moreover, 


dg. 


(o) 


pn 


■dSr 


dv 

r 

dt 


=  V 


df 

r 

ndT 


so  that  Equation  A8  becomes 


and 


g 


(1) 


n 


t   qE    V  df 

v          n  n 

m    "  ,T2  dT~ 

p     V  n 

pn 


pn 


V     g     =  V     [g   (o)   +  eg   (1)    +  0(e2)] 
pnDn         pn  °n  °n  v      ' 


or 


=  v 


t   qE  V     df  m  V  u> 

v^-            n        n  p  p  n  ouqE 

m  V     dt  qE  m 

P            P        n  p 


df 


pn 


V       -   UE    -    T  U)   V      -fir 

n  v  n  n  dt 


which  is  equivalent  to 


dv 

V 

pn 

= 

V 

n 

-  m-E   ■ 

-    T 
V 

n 

The  above  equation  is  consistent  with  steady  flow  results  and  represents 
a  first  order  perturbation  of  the  particle  velocity  due  to  inertia  effects 
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Now,  from  Equation  Ak   the  spectral  drift  velocity  becomes 

df 

vn(uu    ,t)   =  V  f      -V     g     =|iE  +  touV     -J*2- 
DN   n'  n  n         pn°n  v  n  n  at 


The  spectral  slip  then  becomes 


HE 


df 


n 


S        =    £-   +    T    CD       -A- 

n       V  v  vat 

n  n 


But  by  construction,  it  is  reasonable  to  assume  that 


df 

i 

d£" 


ni 


=  1.0 


hence 


where  E  =  EL  is  the  breakdown  field  strength.   Now  V  is  the  amplitude 

at  the  given  frequency  of  the  main  flow  and  it  is  understood  that  where 

V  is  not  significant  (i.e.,  zero)  the  slip  parameter  has  no  meaning. 

Now  p,  and  t   are  a  function  of  the  radius  and  |i  in  addition  of  the 
v 

charge  per  particle  which  in  turn  is  a  function  of  the  radius  and  the 
method  of  charging. 
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Appendix  B 
Effect  of  Particle  Size  Distribution 
The  current  arriving  at  the  probe   is  given  by  Equation  3>   assuming 
one -dimensional  flow. 

00 

Ig(uj)  =  lo  vpKRh(R)~(R)a*  (3) 

where 

00 

f     n(R)dR   =  n 

and  v  is  the  RMS  value  of  the  particle  velocity. 
In  Appendix  A,  it  was  found  that 

dv 

V    =  V   -  U.E  -  T  -T7- 

pn         n  v    dt 

Since  H.E  is   independent   of  time,  under  the  present  assumptions  the 

unsteady  part  of  the  above  equation  is 

f 

i  t  civ 

v       =  v       -  t     — -  (A9) 

pn         n  v    dt 

where  the  prime  denotes  the  fluctuating  components  of  the  velocity. 

Therefore , 

i 

J2E=  1  -  T  i(inv')  (Bl) 

v        v  dt      n 
n 

Now  from  a  linearized  version  of  the  Navier-Stokes  Equation   , 

t 

v  =  V  exp(-iou  t) 
n    n      n  ' 

where 

i  =  imaginary  index 

Thus, 

— (4n  v  )  =  -ioo 
dt v    n       n 
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and  Equation  Bl  becomes 


JS  =  1  +  10)  T 


n  v 
n 

Note  that  the  above  is  time  independent,   so  that  its  RMS  value, 

taking  the  real  part  only,   is 

v  i 

-22  =    (1  -  uu2  t   2)2  (B2) 

v  v  n  v   '  v      y 

n 

Figure  15  shows  a  plot  of  this  function  for  two  arbitrary  frequencies. 

It  can  be  seen  that  because  of  the  strong  dependence  of  t   on  the  radius 

v 

the  function  plots  as  a  low-pass  filter  which  may  be  idealized  for  the 
present  purposes  to  a  step  function. 
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